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ABSTRACT 

In this paper, we show that the eccentricity of a planet on an inclined orbit with 
respect to a disc can be pumped up to high values by the gravitational potential of 
the disc, even when the orbit of the planet crosses the disc plane. This process is an 
extension of the Kozai effect. If the orbit of the planet is well inside the disc inner 
cavity, the process is formally identical to the classical Kozai effect. If the planet's 
orbit crosses the disc but most of the disc mass is beyond the orbit, the eccentricity 
of the planet grows when the initial angle between the orbit and the disc is larger 
than some critical value which may be significantly smaller than the classical value of 
39 degrees. Both the eccentricity and the inclination angle then vary periodically with 
time. When the period of the oscillations of the eccentricity is smaller than the disc 
lifetime, the planet may be left on an eccentric orbit as the disc dissipates. 
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1 INTRODUCTION 

Among the 240 extrasolar planets that have been detected so 
far with a semi-major axis larger than 0.1 astronomical unit 
(au), about 100 have an eccentricity e > 0.3. Five of them 
even have e > 0.8. Such large eccentricities, which cannot 
be the result of disc-planet interaction (Papaloizou et al. 
2001), are probably produced by planet-planet interactions, 
either through scattering or secular perturbation (see Ford 
& Rasio 2008 and references therein), that occur after the 
disc dissipates (Juric & Tremaine 2008, Chatterjee et al. 
2008, Ford & Rasio 2008). 

Here, we show that high eccentricities can be pumped 
by the disc if the orbit of the planet is inclined with respect 
to the disc. The process involved is an extension of the Kozai 
mechanism, in which a planet is perturbed by a distant com- 
panion on an inclined orbit (Kozai 1962). While the Kozai 
effect has always been studied for the case in which the com- 
panion is far away from the planet, the process investigated 
here is shown to be efficient even if the orbit of the planet 
crosses the disc. The classical Kozai effect has of course been 
very well studied. Here, we show that some significant differ- 
ences occur when the classical scenario is extended to apply 
to a disc. 

In section [2] we review the Kozai effect, and show that 
the same behaviour is expected whether the planet is per- 
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turbed by a distant companion or by a ring of material or- 
biting far away. In section [3] we present the results of nu- 
merical simulations of the interaction between a planet on 
an inclined orbit and a disc. We show that, provided most of 
the mass in the disc is beyond the orbit, and the initial in- 
clination is larger than some critical value, the gravitational 
potential from the disc causes the eccentricity and the in- 
clination of the planet's orbit to oscillate with time. This 
may occur even if the orbit crosses the disc. In section [4] we 
summarise our findings, and discuss under which conditions 
this mechanism could operate. The important result is that 
a planet on an inclined orbit with respect to the disc and 
located in or within the planet formation region may have 
its eccentricity pumped up to high values by the interaction 
with the disc. This is of astronomical interest, since inclina- 
tions are beginning to be measured for extrasolar planets. 



2 REVIEW OF THE KOZAI EFFECT AND 
EXTENSION TO A DISC 

We consider a planet of mass M p orbiting around a star of 
mass M* which is itself surrounded by a ring of material 
of mass JV/disc. The ring is in the equatorial plane of the 
star whereas the orbit of the planet is inclined with respect 
to this plane. The motion of the planet is dominated by 
the star, so that its orbit is an ellipse slightly perturbed by 
the gravitational potential of the ring. We study the secular 
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perturbation of the orbit due to the ring. We denote by 
(X, Y, Z) the Cartesian coordinate system centred on the 
star and (r, ip, 6) the associated spherical coordinates. The 
ring is in the (X, y)-plane between the radii Ri and R > 
Ri. We suppose that the angular momentum of the disc is 
large compared to that of the planet's orbit so that the effect 
of the planet on the disc is negligible: the disc does not 
precess and its orientation is invariable. The gravitational 
potential exerted by the ring at the location of the planet 
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where the subscript p refers to the planet and S(r) is the 
mass density in the ring. We assume: 
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with 



wiLii // = Ri/R . We suppose that Ri 3> r p , so that the 
square root in equation (JTJ can be expanded in r p /r and 
integrated to give: 
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In the classical Kozai effect, the planet is perturbed by 
a distant companion of mass M. If we assume the orbit of 
this outer companion is circular of radius R^> r p and lies 
in the (X, y)-plane, then the potential averaged over time 
it exerts at the location (r p ,9 p ) is: 



GM 



1 + 



r 2 P ( 
2i? 2 V 



-1+2 sm ' 



(6) 



Because r p and 9 P appears in exactly the same way in 
$ and <&Kozai, the secular perturbation on the inner planet, 
obtained by averaging over the mean anomaly of its orbit, 
is the same in both cases to within an overall multiplicative 
factor. The results obtained for the classical Kozai effect can 
therefore be extended to the case of the disc. In particular, 
the perturbation due to the disc makes the eccentricity e 
of the planet to oscillate with time if the initial inclination 
angle To between the orbit of the planet and the plane of the 
disc is larger than a critical angle 7 C given by: 
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The maximum value reached by the eccentricity is then (In- 
nanen et al. 1997): 
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and the time t cvo i it takes to reach e max starting from eo is 
(Innanen et al. 1997): 
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with the time r defined as: 

= (l + n)(l-?7-" +2 ) RlNL T 
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where T is the orbital period of the planet and a is its semi- 
major axis. Note that the function (sin 2 7o — 2/5) ~ 1/2 de- 
creases very sharply from infinity to ~ 3 as Jo increases 
from I c to about 45° and then decreases by about 50% as 
7o continues to increase up to 90°. 

The ^-component of the angular momentum of the or- 
bit, L z oc vl — e 2 cos 7, where 7 is the inclination angle 
between the orbit and the plane of the disc, is constant. 
Therefore 7 also oscillates with time and is out of phase 
with e. 



3 NUMERICAL SIMULATIONS 

We consider a star of mass A7* = 1 Mq surrounded by a disc 
of mass Mdisc and a planet of mass M p whose orbit is inclined 
with respect to the disc. The planet interacts gravitationally 
with the star and the disc but we take M p <g Mdi sc so that 
it has no effect on the disc. To study the evolution of the 
system, we use the iV-body code described in Papaloizou & 
Terquem (2001) in which we have added the gravitational 
force exerted by the disc onto the planet. 

The equation of motion for the planet is: 
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where r is the position vector of the planet and $ is the 
gravitational potential of the disc given by equation (fl]) with 
Ri and R being the inner and outer radii of the disc. The 
third term on the right-hand side is the acceleration of the 
coordinate system based on the central star. Tides raised 
by the star in the planet and relativistic effects are included 
through T t ,r, but they are unimportant here as the planet 
does not approach the star closely. Equation (|ll|l is inte- 
grated using the Bulirsch-Stoer method and the integrals 
involved in V$ are calculated with the Romberg method 
(Press et al. 1993). In most runs, the integration conserves 
the total energy of the planet and Lz within 1 to 2%. 

The planet is set on a circular orbit at the distance r p 
from the star. The initial inclination angle of the orbit with 
respect to the disc is 7o. In the simulations reported here we 
have taken n = 1/2 in equation @. The functional form of 
E is shallower than what is usually used for discs, but that 
has no significant effect on the argument we develop here. 

We first compare the numerical results with the analysis 
summarised in section [2] by setting up a case with Ri^> r p . 
In figure Q] we display the evolution of e and 7 for M p = 



10-' M Q , r p 



= 1 au, M d 



= 10" 



M 



0' 



Ro = 100 au, 



Ri = 50 au and 7o = 42°. 3. For this run, Lz is conserved 
within 2% but the energy of the planet is conserved only 
within 10%. We are here in the conditions of the analysis 
of section [2] with r\ — 0.5. From equation Q, we expect 
e ma x = 0.3, which is a bit smaller than the value of 0.41 
found in the simulation. Also the minimum value of 7 should 
be 7 C = 39°. 2 and is observed to be 36°. 5. Note that since 
the energy varies by about 10% in this run, we do not expect 
exact agreement between the numerical and the analytical 
results. We observe that the time it takes to reach e max from 
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the initial conditions is 2.8 x 10 7 years, which agrees well with 
tevoi given by equation ([9]) provided we take e ~ 2xlCT 2 . As 
mentioned above, f evo i becomes very long when 7o is smaller 
than 45°. As the disc lifetime is only a few Myr, e would 
not have time to reach the maximum value in this case, if 
starting from a very small value. 

Figure [2] shows the evolution of e and 7 for M p = 
1(T 3 M , r p = 20 au, M disc = 10 -2 M , R a = 100 au, 
Ri — 1 au and Jo = 47°. 7 (case A). We see that e oscillates 
between e m i n = 10~ 2 and e max = 0.7, whereas 7 oscillates 
between I m - m = 20°. 1 and 7 max = 7o. The values of e m i n , 
e max and 7 m i n differ from those calculated in the analysis in 
section[2] but this is expected as the condition r p <^ Ri, that 
was used in the analysis, is not valid here. However, since 
most of the mass in the disc is in the outer parts, beyond 
the planet's orbit, the behaviour we get here is similar to 
that described in the analysis. The period of the oscillations 
is T osc = 2.2 x 10 5 years. 

For comparison, we have run the classical Kozai case, 
where the disc is replaced by a planet located on a circular 
orbit in the (X, y)-plane. This perturbing planet is at a dis- 
tance R from the central star and has a mass M. We take 
M to be the same as the value of Mdisc above, and consider 
R = 50 and 100 au. The evolution of e and 7 for the inner 
planet in that case is shown in Figure [3] Equation (|8j) gives 
e m ax = 0.5, which is in very good agreement with the values 
of 0.5 and 0.55 obtained from the numerical simulations for 
R = 100 and 50 au, respectively. The time it takes to reach 
e m ax from e m i n , which is T osc /2, is given by equation ((9]) with 
e being replaced by e min and r = [R 3 M*/ (a 3 M)]T / (2tt). 
We get t cvo i = 5.5 x 10 5 and 8 x 10 4 years for R = 100 au 
(e m in = 0.03) and 50 au (e m ; n = 0.02), respectively, which is 
in excellent agreement with the values seen on Figure [3] The 
minimum angle reached when R = 50 au is about 31°, much 
larger than the value obtained when the disc is present. Of 
course this value would decrease if the perturbing planet 
were moved closer to the inner planet, but then the oscilla- 
tions would not be regular anymore, as can already be seen 
in the case R — 50 au. 

Figure [1] shows the evolution of e and 7 for the same 
values of Mdisc and Io as in case A but with M p = 4 x 
10~ 3 Mq, r p — 1 au, R — 50 au and Ri = 0.5 au (case B). 
Here we have e m i n = 5x 10 -2 , e max = 0.72, J m i n = 16°. 4, 
/ mllx = Io and T osc = 2.0 x 10 5 years, comparable to the 
value found in the previous case. 

On dimensional grounds, T osc is expected to be pro- 
portional to r given by equation (|10|l . We have run case A 
with different values of Mdisc and have checked that T osc oc 
1 /Mdisc- The values between which e and 7 oscillate though 
do not depend on Mdi sc , as expected from the analysis. 

We have run case A with different values of r p ranging 
from 2 to 50 au. We observe that for r p roughly below 10 au, 
T sc decreases when r p increases, as expected from the ex- 
pression of t. For larger values of r p though, T OS c does not 
vary much with r p , and is ~ 2-3 x 10 years. For r p = 10 
and 20 au, the extreme values of e and 7 are roughly the 
same if the calculations are started from the same 7o and 
eo, but for r p = 50 au the amplitude of the oscillations is 
very small. In that case, most of the disc mass if not beyond 
the planet's orbit, so that the Kozai effect disappears. 

Finally, we have checked the effect of varying Jo in 
case A. We have found that there is a critical value of 7o, 



7 C ~ 30° , below which eccentricity growth was not observed. 
However, the fact that 7 m i n = 20°. 1 in case A suggests that 
7 C may actually be smaller. Since the time it takes for e to 
grow from very small values when 7q is close to 7 C is very 
long (see equation [9]), the simulations may not have been 
run long enough for a growth to be observed. As 7o increases 
from 7 C to 90°, e max grows from to 1. T OS c is not observed 
to vary significantly with 7o, although the time it takes for 
the eccentricity to grow from its initial value to e max does 
depend on 7o. 

The simulations reported here suggest that when the 
planet's orbit crosses the disc, eccentricity growth occurs 
for significantly smaller initial inclination angles than in the 
classical Kozai effect. 



4 DISCUSSION AND CONCLUSION 

We have shown that, when a planet's orbit is inclined with 
respect to a disc, the gravitational perturbation due to the 
disc results in the eccentricity and the inclination of the orbit 
oscillating if: (i) most of the disc mass is beyond the planet's 
orbit, (ii) the initial inclination angle 7o is larger than some 
critical value 7 C , which may be significantly smaller than in 
the classical Kozai effect. In the simulations we have per- 
formed, in which E oc r _1//2 and Mdisc <C M p (so that the 
effect of the planet on the disc is negligible) , oscillations oc- 
cur as long as the initial planet's distance to the star, r p , 
is smaller than about half the disc radius R . We expect 
a smaller critical distance when E decreases more rapidly 
with radius. Note that 7 C should be independent of the func- 
tional form of E as long as most of the disc mass is beyond 
the planet's orbit. The amplitude of the oscillations depends 
only on 7o. It is small for 7o — 7 C and becomes large when 
7o is increased. In particular, e max — > 1 as 7o — > 90°. The 
oscillations of e and 7 are 180° out of phase. Their period 
Tosc oc 1 /Mdisc When r p < Ro, T osc decreases as r p in- 
creases. For larger values of r p , T osc does not vary much 
with this parameter. In the simulations we have performed, 
Tosc ~ 10 years. As this is much shorter than the disc life- 
time, there is a nonzero probability that the planet is left 
on a highly eccentric orbit as the disc dissipates. 

Note that the process discussed in this paper is not 
a mere trivial extension of the Kozai effect. Indeed, in the 
classical Kozai effect, the periodic behaviour of e is obtained 
because the dependence of the perturbing gravitational po- 
tential on the planet's argument of pericentre cj is through 
a term proportional to cos 2uj. When the orbit of the planet 
crosses the disc, the gravitational potential has a very dif- 
ferent form, and the behaviour of e is much more difficult to 
predict. 

The effect discussed here could be inhibited if other 
processes induced a precession of the planet's orbit on 
timescales shorter than T osc . That may happen if other plan- 
ets are present in the system, or because of dissipative tidal 
torques exerted by the disc, which have been ignored here 
(Lubow & Ogilvie 2001). In the latter case, the Kozai effect 
would then work only for planets orbiting inside the disc 
inner cavity. 

When the planet crosses the disc, loss of energy and 
angular momentum circularises the orbit and aligns it with 
the disc plane on a timescale ~ TM p /(E(r p )R p ), where R p is 
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the planet radius (Syer et al. 1991, Ivanov et al. 1999), which 
can be smaller than T osc for r p larger than a few au. Taking 
this process into account together with the Kozai effect may 
lead to equilibrium values of e and I. This will be studied in 
a forthcoming paper. Note that the usual type II migration 
mechanism that applies to planets orbiting in discs would 
not be relevant here, as it happens when the planet orbits 
in a gap that is locked in the disc evolution. In the coplanar 
case, as the disc spirals toward the central star, it carries 
along the gap and the planet. When the planet is on an 
inclined orbit, it may still open up a gap but it is not locked 
in it, so that it is not being pushed in as the disc spirals in. 

The mechanism described here relies on the planet be- 
ing on an inclined orbit, which could happen as a result of: 

(i) dynamical relaxation of a population of planets formed 
through fragmentation of a protostellar envelope around a 
star surrounded by a disc (Papaloizou & Terquem 2001); 

(ii) mean motion resonances (Thommes & Lissauer 2003, 
also Yu & Tremaine 2001) and (iii) gravitational interac- 
tions between embryos during the planet formation stage 
(Levison et al. 1998, Cresswell & Nelson 2008). In all cases, 
the process that makes the orbits inclined also makes them 
eccentric. According to the results presented here, the disc 
could pump the eccentricities up to even larger values. 

Measures of the projected angle between the axes of the 
planet's orbit and the stellar rotation, using the Rossiter- 
McLaughlin effect, are becoming available. So far, only the 
system XO-3, which has a ~ 12 Jupiter masses planet on a 
3.19 days orbit with e = 0.26, has been shown to have a spin- 
orbit misalignment of at least 37°. 3 (Hebrard et al. 2008, 
Winn et al. 2009). Misalignment has also been reported for 
the system HD 80606, which has a ~ 4 Jupiter mass planet 
on a 111.44 days orbit with e = 0.93 (Moutou et al. 2009, 
Pont et al. 2009). As HD80606 is a component of a binary 
system, the classical Kozai effect could be responsible for 
the misalignment in this system. 

We thank J. Papaloizou and S. Balbus for very useful 
comments which have greatly improved the original version 
of this paper. C.T. thanks the MPIA in Heidelberg, where 
part of this paper was written, for hospitality and support. 
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Figure 1. Eccentricity e (solid line) and inclination angle / (in degrees, dotted line) versus time (in years) for M p = 10 3 Mq, r p = 1 au, 
M disc = 10~ 2 M Q , R = 100 au, Ri = 50 au and J = 42°. 3. 
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Figure 2. Same as figure[T]but for r p = 20 au, Ri = 1 au and Iq = 47°. 7. 
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Figure 3. Classical Kozai effect: Same as figure[2]but with the disc being replaced by a planet of mass M = 10 2 Mq and located at 
distance R = 100 au (upper plot) and 50 au (lower plot) from the star. 
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Figure 4. Same as figure [2] but for M p = 4 X 10 3 Mq, r p = 1 au, R = 50 au, Ri = 0.5 au. 



